Inelastic scattering of light by a cold trapped atom: 
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The light scattered by a cold trapped ion, which is in the stationary state of laser cooling, presents 
features due to the mechanical effects of atom-photon interaction. These features appear as addi- 
tional peaks (sidebands) in the spectrum of resonance fluorescence. Among these sidebands the 
literature has discussed the Stokes and anti-Stokes components, namely the sidebands of the elastic 
peak. In this manuscript we show that the motion also gives rise to sidebands of the inelastic peaks. 
These are not always visible, but, as we show, can be measured in parameter regimes which are 
experimentally accessible. 
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I. INTRODUCTION 



The mechanical effects of photon-atom interactions 
are at the basis of several techniques of manipulation 
of the quantum state of atoms P|. Their understand- 
ing is therefore relevant for implementations of coherent 
control. Features due to the center-of-mass motion have 
been measured in the spectrum of resonance fluorescence 
of trapped atoms |3| and ions 0, 0. IS • These mea- 
surements allowed to gain insight into the state of the 
laser-cooled atoms, and were at the basis of the imple- 
mentation of cooling schemes on the atomic motion based 
on feedback control 8]. 

The spectrum of resonance fluorescence of a laser- 
cooled atom is obtained from the spectral analysis of 
the intensity of the scattered light. By means of these 
scattering processes the atomic center-of-mass motion is 
cooled by transferring center-of-mass energy into the pho- 
tonic reservoir. Hence, the scattered photons carry the 
information about the dynamics that the atom under- 
goes, which can be partly reconstructed in the power 
spectrum. Most studies done so far on the resonance 
fluorescence of trapped ions focused onto the elastically 
scattered radiation by a confined atomic dipole and on its 
Stokes and anti-Stokes sidebands, thereby finding good 
agreement between theory and experiment 2, 3, Q, |a, 0, 
11 H IS E El- On the other hand, the effect of the 
center-of-mass motion on the inelastic part of the spec- 
trum of resonance fiuorescence are largely unexplored. 
Features of mechanical effects in this radiation, which 
has no classical analog, have been studied in 12] for the 
case of a trapped atom, undergoing laser cooling by driv- 
ing two coupled dipole transitions in the regime where 
fluorescence is solely due to the mechanical effects of 
light. In 11 it was found that the quantum motion gives 



rise to sidebands of the inelastic spectrum, which can be 
mapped to Raman processes between the dressed states 
of the bare three-level transition. 

In this work we investigate the spectrum of resonance 
fluorescence of a trapped atom whose dipole transition 
is driven by the cooling laser, and we address in par- 
ticular the issue whether and how features due to the 
mechanical effects of light can be identified and mea- 
sured in the inelastic part of the spectrum. Indeed, we 
show that the quantum motion gives rise to sidebands 
of the peaks of the inelastic spectrum, which have differ- 
ent properties compared to the sidebands of the elastic 
component. These sidebands are not always visible, how- 
ever they can be measured in experimentally accessible 
parameter regimes, like for instance in [l3l |. 

The investigation in this work complements the theo- 
retical study of 10], which focussed onto the Stokes and 
anti-Stokes components of the spectrum of resonance flu- 
orescence of a driven dipole. Like in ^SIj we consider the 
situations when the cooling laser is (i) a running wave 
and (ii) a standing wave. The spectrum is calculated by 
using the spectral decomposition of the atomic master 
equation [ll|ll|ll 

This work is organized as follows. In Scc.|n]we present 
and discuss the spectra of resonance fluorescence. In 
Sec, mil the theoretical model at the basis of the deriva- 
tion is reported. The conclusions are drawn in Sec. IIVI 
The appendices report further details of the derivations 
of Seclml 



SPECTRUM OF RESONANCE 
FLUORESCENCE 



II. 



Let us consider the radiation scattered by a trapped 
atom in the latest stages of laser cooling and detected at 
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FIG. 1: An atom is trapped in a harmonic potential of fre- 
quency V and its dipole transition is driven by a laser field 
with wave vector k and propagating at angle d with the axis 
of motion. The laser cools the atomic motion and a narrow- 
band detector records the spectrum of the intensity of the 
light which is scattered at angle The Gaussian-function de- 
notes the center-of-mass wavepacket, whose finite size affects 
the properties of the scattered light. The inset shows the rel- 
evant internal atomic structure, composed by the ground and 
excited states \g) and |e), which decays at rate 7 into \g). 
Here, Q, and 5 denote the laser Rabi-frequency and detuning, 
respectively. 



of this spectrum. The atomic motion, confined by a har- 
monic potential of frequency couples to the scattered 
light. This gives rise to sidebands of the elastic peak, 
whose origin can be intuitively understood on the basis 
of a classical model for the atomic motion In the 

regime where the atomic wave packet is localized over the 
wave length of light, the mechanical coupling is small 
and scales with the Lamb-Dicke parameter 77 ^ 1. In 
this regime only the Stokes and anti-Stokes sidebands 
are visible. At low temperatures the form of these side- 
bands depends critically on the detection angle. This de- 
pendence is due to interference between quantum paths 
of excitations coupling internal and external degrees of 
freedom |0. 

The inelastic component of the spectrum of resonance 
fluorescence has no classical analog. The mechanical ef- 
fects on this part mainly consists in the appearance of 
sidebands of the inelastic peaks. As we discuss below, 
these sidebands share some properties with but also sub- 
stantial differences from the sidebands of the elastic peak. 



Results 



a narrow-band detector, in the setup sketched in Fig. ^ 
The laser drives the atomic dipole transition, and the 
atomic motion is laser-cooled to a stationary state, which 
is determined by the statistical equilibrium between scat- 
tering processes changing the motional excitations. The 
scattered photons are originated from these processes, 
and the correlation functions and corresponding spectra 
contain some features of these dynamics. The spectrum 
of the intensity of the scattered radiation at the detector 
is 1^ 

/•OO 

S{uj)=xRe dre-''^^(£:(+)(i + r)S(-)(i)) (1) 

where x is a constant, E^^^{t) (E^^^t)) arc the negative 
(positive) frequency part of the electric field at the de- 
tector at time t, and (■) describes the average over the 
atom and electromagnetic field degrees of freedom. For 
suitably chosen detection angles the field at the detector 
is the field scattered by the source, and it is thus deter- 
mined by the state of the atom. 

Let us now make some preliminary considerations, and 
discuss in particular the form of the spectrum when the 
features due to the mechanical effects are negligible. In 
this regime, the spectrum of resonance fiuorescence of 
a dipole transition is constituted by two contributions: 
The elastic component, namely a narrow peak at the fre- 
quency lol of the laser, corresponding to Rayleigh scat- 
tering, and the inelastic component, which originates 
from quantum fluctuations of the atomic dipole. This 
part exhibits the well-known peaks of the MoUow triplet 
which is found when the transition is driven at satura- 
tion [rilfl. 

The mechanical effects of light clearly modify the form 



In this section we report and discuss the theoretical 
results, whose derivation is reported in Sec. lIIII Figure[21 
displays the spectra of the light scattered by the dipole 
transition of an ion which is sideband-cooled (19i. l20j . 
The spectra have been evaluated for different parameter 
regimes. Figures EJa) and (b) correspond to the spec- 
trum at two different detection angles in the situation 
in which the linewidth of the dipole transition 7 — v ji. 
Here, the mechanical effects of atom-photon interaction 
are clearly evident in the sidebands of the elastic peak, 
the two narrow signals at frequency i v. Their width 
scales with the recoil frequency and corresponds with the 
cooling rate 0, |3 0| . These resonances are magnified 
in the inset, where one can see that their centers are 
slightly shifted from the values io j ± v. This shift is due 
to the mechanical coupling [lOl ll^ . The correspond- 
ing curves are a superposition of Lorentz and Fano-like 
profiles, whereby one functional dependence can domi- 
nate over the other depending on the angle of observa- 
tion. The asymmetry of these peaks is due to interfer- 
ence effects, which originate from the quantum nature 
of the ion center-of-mass motion, and whose magnitude 
depends on the detection angle Manifestations of 

the mechanical effects in the inelastic spectrum are small 
but yet visible in Fig. [2tb), causing the appearance of 
an additional resonance on the left part of the curve (at 
Lo ^ uj]^ — 2Av) and a broadening of the central inelastic 
peak. The width of the sidebands of the inelastic peaks 
is substantially the same as the one of the corresponding 
inelastic peak at zero order in the motion. Their height 
scales with the recoil frequency, and, like the sidebands 
of the elastic peak, it is a non-trivial function of the aver- 
age phononic excitation (n) and of the angle of detection. 
Their functional form is reported in Sec. IIII Bl 



3 



■tl} = 40° 




iP = ZOO- 




FIG. 2: Spectrum of the intensity (in arbitrary units) of the light scattered by the dipole transition of a trapped ion in the 
Lamb-Dicke regime. The ion is sideband-cooled by a traveling wave laser. The parameters are 5 — —v, = v, r] = 0.1 and 
(a)-(b) 7 = 0.33!^, corresponding to average phononic occupation (n) = 0.15; (c)-(d) 7 = 0.1;^ ((n) = 0.14). Spectra (a) and (c) 
are measured at the detection angle ip — 40°; Spectra (b) and (d) a.t ip = 200°. The plots show the elastic peak, symbolized 
by the thin line at frequency lu — ujl, the MoUow-spectrum and the motional sidebands of the MoUow peaks. The thin and 
outstanding peaks at about lol i i' are the sidebands of the elastic peak. A magnification of these contributions is shown in 
the insets. 



Figures|2c) and (d) display the spectra of fluorescence 
when the linewidth 7 = v/lO. Here, the mechanical ef- 
fects can be clearly identified in the inelastic spectrum, 
displaced by about ±1/ from the corresponding signal of 
the MoUow triplet. The sidebands of the central inelas- 
tic peak are centered at the same frequency as the Stokes 
and anti-Stokes signals, and give rise to a small broaden- 
ing at their basis. 

We remark that the height of the mechanical sidebands 
- for elastic and inelastic scattering - critically depends 
on the detection angle, as it is visible in the figures. Also 
for the inelastic spectrum this effect is caused by interfer- 
ence of the underlying atomic processes contributing to 
each spectral signature and depending on the direction 
of photon absorption and emission. Thus, the behavior 
also depends on whether the drive is a traveling wave, as 
in Fig. [3 or a standing wave. Figure |3| shows the spec- 
trum of an ion where the trap center is in a point of the 
standing wave (away from the nodes and anti-nodes), for 
7 = 1^/10 and at two different detection angles. Here, 
one can observe the dependence of the sidebands of the 
elastic and inelastic part on the detection angle. Figure01 
displays the spectra of resonance fluorescence at various 
points of the standing wave. For ip = 7r/4, Fig. ^a), 
the spectrum is governed by the outstanding sidebands 
of the elastic peak, overtopping the MoUow triplet. At 



(p = Btt/S, Fig. ^b), the trap is closer to the node of 
the standing wave. Here, the resonances of the MoUow 
triplet get closer. In this case, the height of the side- 
bands of the elastic peak decrease and the sidebands of 
inelastically scattered light are visible. A peculiar case is 
found when the atom is placed at the node. Fig. ^c) at 
ip = tt/2. Here, the only spectral signals are the motional 
sidebands. These sidebands are now perfect Lorentzian 
curves of identical shape independent of the detector po- 
sition, as already pointed out in What is novel is 
that all other contributions of the spectrum, elastic and 
inelastic, disappear due to the vanishing of the driving 
field at the center of the trap. In particular, even the 
motional sidebands of the inelastic spectrum disappear. 
The appearance of the sidebands of the elastic peak is 
due to the effects of the gradient of the field intensity 
over the atomic wavepacket at the node. The suppres- 
sion of inelastic scattering processes, included the ones 
that change the motional excitation, can be understood 
by considering that there are no quantum fluctuations of 
the atomic dipole at the node of a standing wave. 

These behaviors are discussed in detail in the follow- 
ing theoretical treatment. It should be remarked that 
the sidebands of the inelastic component become visi- 
ble when the transition, at zero order in the mechani- 
cal effects, is driven at saturation. Below saturation the 
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FIG. 3: Spectrum of the intensity (in arbitrary units) of the 
light scattered by the dipole transition of a trapped ion in the 
Lamb-Dicke regime. The ion is sideband-cooled by a standing 
wave laser. The parameters are 5 — —v, f2 = !/, 7 = O.li/, 
77 = 0.1 and detector angle (a) i/- = 40°; (b) i\) = 200°. The 
standing wave is described by the function cos(fca; + kp) where 
(/3 = 7r/4 is the value at the trap center. The mean phononic 
excitation is (n) = 0.04. The insets display the sidebands of 
the elastic peak, the dashed-dotted lines shows the inelastic 
spectrum contribution. 
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mechanical effects in the inelastic spectrum manifest in 
small sidebands of the Lorentz curve centered at the laser 
frequency. These sidebands are at the same frequencies 
as the ones of the Stokes and anti-Stokes components of 
the elastic peak, and are covered by these signals. 



FIG. 4: Same as Fig.j^for different position in the standing 
wave. Here, (a) -p = 7r/4 ((n) = 0.04); (b) ip = Svr/S ((n) = 
0.002); (c) = 7r/2 ((n) = 0.0006). The parameters are 
5 = -z/, n = 7 = O.ljy, r] = 0.05, V = 200°. The insets 
display the sidebands of the elastic peak, the dashed-dotted 
lines correspond to the inelastic spectrum contribution. 



III. EVALUATION OF THE SPECTRUM OF 
RESONANCE FLUORESCENCE 

We consider an atom of mass M, whose center of mass 
motion is confined by a harmonic potential of frequency 
V. We assume that the relevant motion is along one di- 
mension, the X-axis, while the radial motion is frozen out. 
A laser, of frequency and wave vector fc, couples quasi- 
resonantly with the atomic dipole transition between the 
ground and excited states \g) and |e) at frequency and 
linewidth 7. The motion is in the Lamb-Dicke regime, 
namely the size of the center-of-mass wave packet is as- 
sumed to be much smaller than the laser wavelength, and 
the dipole linewidth ^ < v. In this regime the laser can 



be set at freque ncy ojl — ujo — v thereby sideband-cooling 
the motion [l9l |2(]| . The atomic motion is in the station- 
ary state of sideband cooling and scatters the laser light. 

Let a detector measure the radiation scattered at the 
angle i/; with the atom motional axis. In the far-field the 
scattered field is proportional to the retarded value of 
the atomic dipole at the corresponding position x. The 
spectrum at frequency uj is given by |2lj| 

S{u;) = xWRe / dre-'^"^-^^)^ {D^ (x, t)D{x, 0)U 
Jo 

(2) 

where x(V') is a function of the emission angle and D{x, t) 
is dipole lowering operator in the reference frame rotat- 
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ing with the laser frequency ojl at the position x of the 
atomic center of mass. The correlation function in Eq. (O 
is formally evaluated by means of the quantum regression 
theorem, according to which D{x,t) = D{x)e^'^ with C 
the Liouvillian of the master equation for the atomic den- 
sity matrix g, 

dg/dt = Cg. (3) 

The average (•)st is taken over the atomic density matrix 
Pst at steady state, solution of the equation Cgst = 0. 

The expHcit form of S{u}), Eq. jSJ, is found by applying 
the spectral decomposition of the Liouville operator C 
according to the secular equations 

Cg^ = \g^, 
g^C ^ Xg^, 

with eigenvalues A and right and left eigenelements p'^ 
and g^, respectively [Ij, [13 • The orthogonahty and 
completeness of the eigenelements is defined with respect 
to the trace, such that TT{g^' g^} = S^''^. We define 
the projectors onto the eigenspace corresponding to the 
eigenvalue A as — g^ (i^ g^ such that their action on 
an operator X is defined as V'^^X = g^Tilg^X}, and 
they thus fulfill CP^ = V^C = XV^. By applying this 
formalism, we rewrite Eq. (0) as 

S{lu) = x(V')Re V ^^^Tr {D^x)V^Dix)g,,} . 

(4) 

where we have used the completeness relation of the 
eigenelements of C. Although completeness of this kind 
of basis is generally not warranted, the spectral decom- 
position that we consider in the following is complete, 
and the spectrum of resonance fluorescence can be thus 
cast in the form of Eq. |0J . 

A. Model 

We introduce now the basic elements determining the 
dynamics of the atom. The coupling with radiation is 
assumed to be in the Lamb-Dicke regime. This regime 
is characterized by the small value of the Lamb-Dicke 
parameter 

which scales the coupling between internal and external 
atomic degrees of freedom due to photon scattering. We 
identify 77 as the parameter of the perturbative expan- 
sion. According to this expansion, at second order the 
spectrum in Eq. Q is decomposed into the terms 

5H = S'oH + 5iH + 52H + 0(?7^) (6) 

where the subscripts label the corresponding order in the 
Lamb-Dicke expansion. In order to evaluate 5*^(0;), we 



consider the operators C and D at second order in the 
expansion in 7y. The dipole = Do + Di + D2 + 0{r]'^), 
where the individual terms have the form 

Do=(J 

Di — — ifcx cos '0 cr, 
D2 = -ifc^x^ cos^ipa, 

with a ~ \g){s\- The Liouville operator is decomposed 
into the terms C = Co + Ci+ C2, which we now introduce 
in detail. At zero order internal and external degrees of 
freedom are decoupled, namely 

£0 = A + (7) 

where £e and Ci act on the external and internal degrees 
of freedom, respectively. Here, 

-CeP = r^[i?moc, £»] (8) 

in 

with 

iJmcc = fti' (^a^a+ , (9) 

where a and a' are the annihilation, creation operators 
of a quantum of energy hu, respectively. We denote 
with |n) the eigenvectors of ffmec, fulfilling Hynedn) = 
hv{n + l/2)|n) with n = 0, 1, 2, . . . The term L\ acts on 
the internal degrees of freedom and it is defined as 

Cig^hH^ + Va,g\+lCg (10) 
in 

where i/g = f)-5\g){g\, with 5 = lul — t^o, and 

K-og — (2aga^ — cr^ag ~ ga^cr) 

are the interaction with radiation at zero order in the 
mechanical effects. Here, we have denoted by ^Ci'P) the 
laser Rabi frequency, whereby C{ip) is a dimensionless 
function of the angle ip, which is Cif) = ^^p{^f)\v=o for 
a traveling wave drive, and is given by Cif) = cos for 
a standing wave drive. 

The first and second order Liouvillian give rise to the 
mechanical coupling between internal and external de- 
grees of freedom, 

Cig = ^[xVi,g], (11) 

C2g - ^ [2:^2, e] +JC2g. (12) 

where 

Vi = ^kcosenC{(p)<T'' +R.C. (13) 

V2 = ifc2 cos^ enc'i'py'' + h.c. (m) 

f^2P — P^k'^cr (2xpx — x^p — px^) (15) 
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Here, 9 is the angle between the laser and the trap axis, 
/3 is a constant, which gives the average recoil due to 
spontaneous emission of a photon, and C'(¥')i C"(¥') are 
the first and second derivative of C(v) with respect to ^p. 

The spectral decomposition of the Liouville operator is 
obtained using the eigenvalues and eigenelements of op- 
erator £o : of which we evaluate the corrections at second 
order in "q. We denote the eigenvalues and eigenelements 
of Cq by Ao and g^, Qq, respectively. Since internal and 
external degrees of freedom are decoupled at zero order 
in T], the eigenvalues of Cq are Aq = Ai + Ae, and the 
eigenelements 



Qo 



where Cip^^ — Xip^^ and -CeA*'^'^ = Ae/^"^^- Correspond- 
ingly, the projector into the subspace at the eigenvalue 
Ao is 

whereby their action on the operator X is defined as 

r^'X = p^'TtiIp^'X} 
U^'^X = //>=TrE X} 

and Tri (Tte) denotes the trace over the internal (exter- 
nal) degrees of freedom. 

The spectrum of Ci characterizes the dynamics of the 
two-level transition and the spectral properties of the ra- 
diation emitted by the bare atom The eigenvalues 
of £e take on the values Ae = itv, with £ = 0, ±1, ±2, . . . 
Each cigenspace at Ae is infinitely degenerate, and the 
corresponding left and right eigenelements are, for in- 
stance, jil^ = \n + t){n\, /i^ = \n){n + £\. These eigenele- 
ments constitute a complete and orthonormal basis over 
the eigenspace at this eigenvalue. In particular, the pro- 
jector over the eigenspace at Ae = iii^ is defined on an 
operator X as 

U^^-^'^X = ^/„TrE{/ilX} (16) 

n 

= ^\n){n\X\n + £){n + i\, 

n 

where Trs denotes the trace over the external degrees of 
freedom. 

At higher orders in the expansion in 77, internal and 
external degrees of freedom are coupled, and the degen- 
eracy of the subspaces at eigenvalue Ae is lifted 0, |22| . 
The perturbative corrections to the eigenvalues Aq, to the 
eigenelements pg , Pq, and to the projectors Vq are found 
by solving iteratively the secular equations at the same 
order in the perturbative expansion, and are given by 



where Vn is the zero-order projector onto the subspace 



at eigenvalue A, 'Pq = ® ^o- Using ((T7|l we obtain 



(1 



K 

(Ai 



^92 = 



A) 



An — Cf) 



(18) 



Ao — jOq 



(Ai - A) + (A2 - £2) 



£'0- 



Analogously, one finds the perturbative corrections to the 
left eigenelements . This in turn allows one to evaluate 
the perturbative corrections to the projectors Vq, which 
are immediately found by using the explicit form of Qj, 
into the relations 



00 92 



A -A 

9i9o 
9i9i 



92 90 



(19) 
(20) 



The equations for the corrections Ai, A2 to Ao are 

Ai - Tr{g^Aeo} = 0- (21 
A2 = TT{g>;C2P^} + Tr{p>;/:ip^}, 

= Tr{p^C2P^} + TT{p^Ci 



^—^C^p^}, (22) 
^0 ~ i-O 



where we have used relation V^CiV^ = in Eq. (EU- 

A relevant eigenelement of this spectral decomposition 
is the one at Ao — 0, namely the steady-state density 



matrix. At zero order this is given by Qq = poP-, where 
p is the density matrix for the external degrees of freedom 
in the final stage of the laser-cooling dynamics and has 
the form 



1 



in) 



l + {n) 



(23) 



and 



= Tr{a^a^} 



(24) 

is the average phonon number at steady state. Corre- 
spondingly, ^g^" = IiIe, where Ij is the identity for 
the Hilbert space of the internal (j = I) and the external 
{j — E) degrees of freedom. 



B. Explicit form of the sidebands of the elastic and 
inelastic components 

We now proceed in evaluating the spectrum of res- 
onance fluorescence using the terms introduced in the 
previous section. At zero order in the mechanical effects, 
the spectrum is 



1 



(lu - ujl) - A 



G(A) 



(25) 



i^-r^)9'i 



Ao — Cq 



(Ai -£l)^^^, 



(17) 



where G(A) = TtID^Vq DqPq} , which corresponds to 
the spectrum of the bare atomic dipole. It thus exhibit 
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peaks which are located at frequencies equal to the real 
part of the eigenvalues Aj, whose width is determined by 
the imaginary part of Aj. Note that TrlZ^JPg Z?oPo} = 
TT{a^ PQ^}TT{pQ^apo}. For later convenience, we intro- 
duce the scalar functions 

g(Ai) = Tr{/5^Vpo}, 

such that with this notation G(A) — g{\i)g{Xi). 

The term S'i(w) = 0, since the motional steady state 
p, Eq. if^ . is a thermal state, and thus diagonal in the 
vibrational number basis |l2j . 

The features due to the mechanical effects thus man- 
ifest at second order in the Lamb-Dicke expansion, and 
the corresponding term has the form 



{u - wl) - A 



F(A) 



where 



F(A) 



^ Tt{dIv^d,p,} 



(26) 



(27) 



a+b+c+d=2 



The function F{\) can be decomposed into three contri- 
butions, namely 



F(Ae,Ai) = Fo(Ai)5ak,o 

+ F+(Ai)(5ae,,, + F_(Ai)5ab,- 



(28) 



This decomposition allows one to identify the effects of 
the quantum motion on the features of the spectrum. 
The term Fo(Ai) is a second-order correction to the fea- 
tures of the zero order spectrum, So{oj). Namely, it gives 
rise to small corrections to the curves of the zero-oder 
spectrum, however it does not affect relevantly its gen- 
eral form. The terms -F±(Ai), on the other hand, give rise 
to novel peaks centered at the frequencies Im{Ai} ± v, 
which can be identified with the sidebands of the peaks 
appearing in Sq{lu). These terms have the form 



and 



F+(Ai) —rj'^ cos^ 6 r(Ai, zi/) {r*(Ai, w) 4- m(Ai, w)} {n) - t(Ai, iu) 



-rj cos 6* cos V' 



g(Ai)[{r,(Ai,ij/) + u{Xi,w)} (n) -t{Xi,ii')] + r{Xi,ii')g{Xi){n) 



+Tf cos^ g{Xi)g{Xi){n) 



F-{Xi) =Tfcos^e r{Xi,~iv) {r^{Xi, -iv) + u{Xi, ~iv)} {{n) + l)+t{Xi,-iv) 



-rf cos 6 cos ?/> 



.g(Ai) [ {r,(Ai, -iv) + ?i(Ai, --iv)} {{n) + 1) + t(Ai, -iv)] + r(Ai, -iv)g{Xi){{n) + 1) 



-ff cos2V<7(Ai)5(Ai)((n} + l), 



(29) 



(30) 



where we have used rj = kxo with xo — \/ fijlmv and we 
have introduced the functions 

r(Ai, Ae) = iTri{i?5(Ai + Ae - A)"Ml4, Po^]} (31) 

r.iXi.M) = iTri{p^'i?o(AE + /:i)"M^i,Po]} (32) 

u(Ai, Ae) = -■^Tri{p^^[yi, (Ai + Ae - A)~'^oPo]} 

(33) 

i(Ai,AE) = iTri{p^i?o(AE + A)"Vo^^i} 

- ■^Tri{Fip^(AE + Ai - UY^D^p^}. (34) 

The basic steps that lead to Eqs. (|29|I - H34|) are reported 
in App.^ The results are plotted and discussed for some 



parameter regimes in Sec. Ull 

We now discuss some general properties of the equa- 
tions for two cases: When the laser, driving the dipole 
and cooling the motion, is a traveling wave and when it 
is a standing wave. 

Traveling wave drive. The case when the driving laser 
is a traveling wave is found by setting C'iv) = i at (/? = 
in Eq. We evaluate the terms Eqs. (jSIJ - lESI) using 
the zero-order eigenelements p^ and pp determined by 
the eigenvalue equations CoPi — ^oPi and PiCo — Aop^, 
as in [l^. The eigenvalue Ai = gives the contribution 
of the elastic peak. As we focus on the case of a dipole 
driven at saturation, the three other - nonvanishing - 
eigenvalues of the zero-order spectrum determine the po- 
sition and the widths of the MoUow peaks. 

We characterize the higher-order contributions to the 
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spectrum by their corresponding eigenvalues at zero order 
and analyze the single terms composing 5*2 (tf), Eq. (|26|l . 
As 5*0(0;) never vanishes for traveling wave drives, the 
term Fq gives rise to small corrections to the elastic peak 
(component at Ai — 0) and to the MoUow spectrum (com- 
ponents at Ai 7^ 0). The terms ^±(A) correspond to novel 
peaks which add to the ones of the zero-order spectrum. 
At Ai = these terms give rise to the motional side- 
bands of the elastic peak, centered about Ae = ±ii/. In 
this case the terms Eqs. (IST|) - have to be evalu- 
ated for Pq' = li and p^^ ~ pQ. By applying the cyclic 
properties of the trace one finds that u(0, Ae) = and 
that r*(0, Ae) — r'(0,AE)*- Since Ai = 0, the width and 
shape of these curves is critically determined by the sec- 
ond order correction A2 to the eig envalue. The explicit 
evaluation can be found in |l2ll22j | The result we find is 
in agreement with the one reported in [Toj . 

The mechanical effects on the inelastic part of the spec- 
trum manifest themselves in the components at Ai 7^ 0. 
They give rise to small shifts of the center-frequency 
Im{Ai} and of the linewidth Rc{Ai}, and to novel peaks 
which are centered at the frequencies ImAi ± v, and 
which we discuss in the following. These peaks are 
the sidebands of the corresponding zero-order curves in 
the MoUow spectrum. Their width is the width of the 
corresponding zero-oder MoUow peak, namely Re{Ai}. 
The intensity of the particular signals is proportional to 
F±{Xi), where generally all the terms - g{Xi) 
and g{Xi) contribute. Their shape is in general a su- 
perposition of Lorentz (real part of F±(Ai)) and Fano- 
like (imaginary part of i^±(Ai)) curves. Note that in the 
plots shown in Sec.m the real part of F±{Xi) dominates 
for the visible peaks and thus these sidebands appear as 
Lorentz curves in the spectrum. Their visibility, how- 
ever, depends on the detection angle. Indeed, also for 
these terms there are interference effects that give rise to 
an asymmetry between the right and the left sideband 
of a given inelastic peak. Note that the sidebands of the 
central inelastic peak are located at the same frequencies 
of the Stokes and anti-Stokes sidebands. These terms 
do not vanish, but are usually hidden by the Stokes and 
anti-Stokes contributions. Spectra for different parame- 
ter regimes are shown in Fig. [3 

Standing wave drive. The case when the driving laser 
is a standing wave is found by setting ((ip) = cos(p in Vq 
and Ci'P) = ~ sin 95 in Eq. H13|l . The general characteris- 
tic of the terms composing the spectrum is basically the 
same as for a traveling wave, the difference arises mainly 
in the dependence of the curves on the scattering angle. 
The analytical form of the two sidebands, which we ob- 
tain by evaluating explicitly i^±(0) and the linewidth A2 
of the resonances, is in agreement with the formulas re- 
ported in [lOi] . A peculiar behavior is found if the atom 
is placed at the node of the standing wave, ({(p) — 0, 
which we discuss in the following. Indeed, in the node of 
a standing wave the electromagnetic field vanishes at zero 
order in the Lamb-Dicke expansion. Hence, the dipole is 
excited by processes at higher order in the Lamb-Dicke 



parameter 77. In other words, in this regime the spec- 
trum of resonance fluorescence is only due to the gradi- 
ent of the electromagnetic field over the finite size of the 
ion wave packet, and thus originates from mechanical ef- 
fects of atom-photon interaction. In the formulas the 
terms Vq — V2 — and the steady-state density matrix 
Po = 1.9) (sIm- Details over the eigenelements and eigen- 
values of Ci are reported in App.0 In this case DqPq = 
and consequently the term u{X) = 0, see Eqs. (|33|) . For 
the same reason the second term on the RHS of t{X) 
vanishes, see Eq. H34|l . Also the first term on the RHS 
of t{X) vanishes, as one can verify by taking the explicit 
form of Vi, Eq. (|13|l . and of the eigenelements. Hence 
t(A) = 0. Moreover, in Add. iBlwe show that Fq{X) — 0. 
The only terms which do not vanish are r(0,AE) and 
r40,AE) in F±(A), Eqs. Therefore, the spec- 

trum is at second order in the Lamb-Dicke expansion, 
and the only non- vanishing contributions are the Stokes 
and anti-Stokes sidebands of the elastic peak. Spectra 
evaluated for different parameters are shown in Figs. El 
and El 



IV. CONCLUSIONS 

The spectrum of resonance fluorescence of a cold 
trapped ion, which is in the stationary state of laser cool- 
ing, presents several features due to the mechanical ef- 
fects of light on the quantum motion of the atomic center 
of mass. In this work we have shown that, in addition to 
the well-known Stokes and anti-Stokes components, sim- 
ilar features appear also in the inelastic spectrum. These 
features are sidebands of the inelastic peaks and have no 
classical analog. The quantum features in the mechani- 
cal effects of atom-photon interaction become visible at 
low temperatures, when the motional state is close to 
the ground state, and when the transition linewidth is 
smaller than the trap frequency, namely in the regime 
where the ion can be sideband cooled. This parameter 
regime can be encountered, for instance, in experimental 
situations like 

In our investigation we have considered the light scat- 
tered from a laser which sideband-cools the motion in 
two cases: when the drive is a traveling wave and when 
it is a standing wave. A peculiar behavior is observed 
when the drive is a standing wave: Here, depending on 
the position of the ion in the mode, the spectrum changes 
dramatically. In particular, in the node of the standing 
wave we find that the spectrum is solely composed by the 
Stokes and anti-Stokes sidebands, while any other contri- 
bution is at higher order in the Lamb-Dicke expansion. In 
this regime we recover the form of the sidebands compo- 
nents as predicted in |lClj |. The situation in the node can 
be compared with a A-system driven at two-photon res- 
onance, where the zero-order contribution vanishes due 
to the existence of the dark state coherence. In that 
case, however, sidebands and unshifted contributions of 
the inelastic spectrum are visible, since the electric field 
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does not vanish at zero order in the Lamb-Dicke expan- 
sion [13. 

In conclusions, we have studied the effects of the quan- 
tum center-of-mass motion on the inelastic part of the 
spectrum. This work complements previous theoretical 
analysis which investigated the elastic part. In 

the future we will look at how the quantum center-of- 
mass motion affects higher order correlation functions of 
the scattered light. 
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APPENDIX A: EVALUATION OF F±(A) 

The terms of Eq. H28() which contribute to ^±(A) are 

/i(A) = Ti{Dlr^Dogi} + Tr{Dlv^Dogo} 
/2(A) = Tv{Dlv^Dogi} + Tr{Dlv^Dogo} 
/3(A) = Tr{Dlr^D,go} 
/4(A) = Tr{D\r^D,go} 

We evaluate now their explicit form using the perturba- 
tive expansion using the definitions introduced in Sec. lIIII 
We denote by the superscript * the part of the term, 
which contributes to F±{X). The relevant part of the 
first term can be rewritten as 

, cos^ e 

J I (A) = r(Ai,AE) 

X [TrE {a; U""^ a;//} (Tn {p^' Z?o ( Ae + A ) " ' [14 , Po] } 

- Tri{p^' [Vi, (Ai + Ae - Cir^DoPo]}) 

- Tri{yip^^(Ai + Ae - A)"'^oPo}^ 



obtain 

{r,(Ai, Ae) + m(Ai, Ae)} Ti^ixU'^^xp.} 
-f i(Ai,AE) Ty^{xU^^[x,p\} . 

It can be simply verified that function r*(Ai,AE) in 
Eq. fulfills the relation r^{{),X^) = r(0, Ae)*. More- 
over, the term u(0, Ae) = 0. Finally, the function 
t(Ai,AE), defined in Eq. 1^1)1 . reduces for Ai = to the 
function t{\^) in [1(tI |. 

Using these definitions, the other terms can be rewrit- 
ten as 

/^(A) = cos6'cos7/' 5(Ai)x 

{n(Ai, Ae) + u(Ai, \^)}l:Y^{xU^''x^l} 
^t{Xi,X^)Tr^{xU^^[x,p\} 

and 

/3±(A) = -k^coaOcos^l^ 5(Ai)r(Ai, Ae) Tv^ixU^'^x^i} 
ft[X) = fc^cos^V' <7(Ai)5(Ai) Ty^ixU^^xtx} 

where we have used the definitions from Eq. (|26|l . The 
trace terms over the external degrees of freedom are con- 
veniently evaluated using the basis set corresponding to 
the projectors in Eq. H16|l . giving 



Tr^ixU'-^xp.} = xl [5x^,-^,{{n) + 1) + 6x^,,,{n)] 
Tve,{xU^'^[x,^l]} = xl [-~5xe,%u + 



(AI) 
(A2) 



with xq = ^JhjlMv. Using these results into Eq. H28(l 
we obtain Eqs. and 1(201). 



APPENDIX B: ION IN THE NODE OF THE 
STANDING WAVE 

In this appendix we evaluate explicitly the contribution 
-Fb(Ai) and show that it identically vanish in the case 
where the atom is at the node of the standing wave. 

Let us first discuss the terms contributing to Fo(Ai). 
In the node of the standing wave po — I 5 ) ( 5 I ■ By using 
the relation DaPn = (a = 0,1,2 =), the terms which 
do not trivially vanish are 



fo(A/) = /°(Ai) + /5"(Ai) + /°(Ai) 
where are the components of the terms 



where we have used the function r(Ai,AE) defined in 
Eq. H31|) . Note that for Aj = it reduces to the function 
-r(AE) in llo|. Using definitions W^ . and we 



/i(Ai) 
/5(Ai) 

/6(Al) 



T^{DlviD^g^} 
Tr{DjPo^I?oe2} 
TT{DlvlD^g^} 



(Bl) 



(B2) 
(B3) 
(B4) 
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which contribute to Fo(Ai). Their expHcit form is found 
with the internal eigenelements, which at the node take 
the form 



5)(.9l, 
e)(e| 

P' = |e)(5l, 
= l5)(e|, 



'){9\ 



P^* = 1e, 
P°-|e)(e|, 
P+ = l<?)(e| 
p-^\e){9\ 



Using the expressions for the external traces, Eqs. (|A1|) 
and ljA2(l . and the explicit form of eigenelements and 
eigenvectors, one finds that /i(Ai) = — /"(Aj). Thus, 
the two terms mutually cancel in Eq. ljBl|) and hence 
Fo(Ai) = 0. 



with the corresponding eigenvalues A*^' = 0, Aj = —7 
and \i = —7/2 ± i5. Substituting into Eq. HB4|) one 
finds /e (A) = 0. The other two terms give at Ae = 



fc2 COs2 e 



Ay(A+ + Ay 



It{^i) = hi,x- X 



E 



\TI^{xU^'^x^l} i:Yji{xU^'^ ^lx} 



L A"(A+ + Ay AO(A- + A^) 
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